Thermodynamics Predicts How Confinement Modifies Hard-Sphere Dynamics 
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We study how confining the equihbrium hard-sphere fluid to restrictive one- and two-dimensional 
channels with smooth interacting walls modifies its structure, dynamics, and entropy using molecular 
dynamics and transition-matrix Monte Carlo simulations. Although confinement strongly affects 
local structuring, the relationships between self-diffusivity, excess entropy, and average fluid density 
are, to an excellent approximation, independent of channel width or particle-wall interactions. Thus, 
thermodynamics can be used to predict how confinement impacts dynamics. 
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Tlie molecular dynamics of fluids confined to small 
spaces can differ significantly from the bulk. These dif- 
ferences have generated wide interest because confined 
fluids feature prominantly in both nature and technol- 
ogy. Examples include dynamics of water near proteins 
or in concentrated cellular environments, transport pro- 
cesses across biological membranes, and fluid flows en- 
countered in micro- or nanofluidic devices, to mention a 
few. Given that a significant fraction of the molecules 
in these systems populate highly inhomogeneous inter- 
facial environments, it is easy to appreciate why con- 
finement has nontrivial consequences for their transport 
coefficients (e.g., diffusivity and viscosity). Nonetheless, 
a theoretical framework that can reliably predict these 
consequences has been slow to develop. 

One logical starting point is to ask whether confine- 
ment induced modifications to equilibrium fluid proper- 
ties, such as the density, can explain some of the observed 
differences in dynamics fl, 2, 3]. For instance, if the pres- 
ence of a strongly attractive substrate increases the local 
fluid density near the fluid-substrate interface, one might 
naturally expect a corresponding decrease in particle mo- 
bility near that interface, and vice versa. This type of 
argument is physically intuitive, and it has been recently 
used to rationalize why nanoscale materials exhibit glass 
transition temperatures that are shifted relative to their 
bulk values [Sj. However, it seems doubtful that aver- 
age structural quantities alone can account for the wide 
variety of dynamical behaviors observed in both simu- 
lations and experiments of confined fluids Him. As 

a result, it is natural to ask whether other equilibrium 
measures, such as the entropy, can provide additional in- 
sights. Unfortunately, it is currently difficult to obtain 
the necessary experimental data for testing these possible 
connections between thermodynamics and dynamics for 
confined fluids. Thus, simulation results on simple and 
well-defined model systems are of great complementary 
value. 

In this Letter, we advance the current understanding 



of the relation between thermodynamics and dynamics 
in inhomogeneous systems by resolving the following two 
questions, (i) Can either the entropy or the average den- 
sity be used to determine the extent to which confinement 
alters the diffusivity of a hard-sphere (HS) fluid? (ii) 
If so, do the specific interactions between particles and 
the channel boundaries significantly impact the result? 
While the confined HS system represents arguably the 
most elementary and well-studied model for inhomoge- 
neous colloidal and molecular fluids, there is still surpris- 
ingly little known about the possible connections between 
its basic thermodynamic and kinetic properties. If such 
connections do exist and prove to be robust, it suggests 
that equilibrium theories for inhomogeneous fluids might 
generally provide important information regarding how 
confinement modifies the transport properties of fluids. 

One reason to speculate that entropy could be a re- 
liable predictor for how confinement affects the diffu- 
sivity is its empirical success for capturing the dynam- 
ical behavior of bulk materials. In particular, computer 
simulation studies have demonstrated that the single- 
component HS and Lennard- Jones fluids, along with a 
variety of models for liquid metals, exhibit, to a very 
good approximation, a one-to-one relationship between 
diffusivity and excess entropy over a broad range of ther- 
modynamic conditions 0j H S 0| • Excess entropy 
has also been shown to accurately capture the behavior 
of diffusion phenomena in ffuid mixtures 0, lisl as 
well as those in solid-state ionic conductor and quasicrys- 
talline materials |^]. 

To explain the origin of the correspondence between 
excess entropy and diffusivity in bulk materials, several 
researchers have presented independent derivations of ap- 
parent scaling laws relating the two quantities. The ear- 
liest that we are aware of is due to Rosenfeld and is 
motivated by a variational thermodynamic perturbation 
theory p, |2| . Dzugutov later used arguments based on 
kinetic theory to justify a similar scaling [sj, and recently 
mode-coupling theory has been employed to establish an 
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approximate basis for the observed connection 0, ■ 
Despite the large amount of effort that has been devoted 
to justifying these scaling laws theoretically and testing 
their validity for bulk materials, to our knowledge the 
relationship between excess entropy and diffusivity has 
never been tested in inhomogeneous fluids, nor has it 
been used as a tool to understand how confinement af- 
fects dynamics. 

To carry out such a test, we studied how the structure, 
thermodynamics, and dynamics of the single-component 
HS fluid confined to restrictive two-dimensional {2d) or 
one-dimensional (Id) channels ("films" or "pores", re- 
spectively) bounded by smooth, interacting walls differ 
from those of the bulk system. We considered five dif- 
ferent 2d channel sizes that were effectively macroscopic 
in the x and y directions and had particle-center acces- 
sible dimensions in the confining z direction of hz — 15, 
10, 7.5, 5, and 2.5, repectively. We also considered three 
Id channel sizes that were effectively macroscopic in the 
X direction and had particle-center accessible dimensions 
in the confining y and z directions of hy x = 7.5 x 7.5, 
7.5 X 5, and 5x5, respectively. To simplify the nota- 
tion, we have implicitly non-dimensionalized all lengths 
in this study by the HS particle diameter a and all times 
by the combination a^/mfi, where m is particle mass, 
(3 — l/fceT, fee is Boltzmann's constant, and T is tem- 
perature. Consequently, all energies are given per unit 
k^T . Position-dependent interations between the parti- 
cles and the confining channel walls u{C,) were calculated 
using a square-well potential: 
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where C, represents the shortest distance between a given 
particle center and the wall of interest, and Cw is the 
strength of the effective particle-wall interaction. We 
considered five specific particle-wall interactions for the 
2d channels: e„ = 1 and ew = 0.5 representing repulsive 
walls, Cw = representing "hard" but neutral walls, and 
Ew = —0.5 and e„ = — 1 representing attractive walls. 
Only the hard walls were considered for the three highly 
restrictive Id channels. 

To monitor kinetic processes in these various systems, 
we performed a series of event-driven molecular dynam- 
ics simulations jl5j in the microcanonical ensemble using 
N — 4500 hard spheres. Periodic boundary conditions 
were employed in the d "free" directions (i.e., directions 
not confined by walls). The dimensions of the simula- 
tion cell in the periodic directions were set to various 
values to simulate fluids with different average number 
densities that span the stable equilibrium range, from 
the dilute gas to the fluid at its freezing transition. We 
extracted the self-diffusivity of the fluid D by fitting the 
long-time (i 3> 1) behavior of the average mean-squared 



displacement of the particles to the Einstein relation 
< Ar^ >= 2dDt, where Ar^ corresponds to the mean- 
square displacement in the d periodic directions. We also 
calculated D for several state points with both smaller 
{N = 3000) and larger {N = 6000) particle numbers to 
verify that system size effects in the periodic directions 
were negligible. 

We determined the behavior of the excess entropy 
per particle s'^^ using grand canonical transition-matrix 
Monte Carlo (GC-TMMC) simulations [3. Here, s'''^ 
is defined to be the difference between the entropy per 
particle of the fluid and that of an ideal gas with the 
same spatial distribution of the particle number den- 
sity. GC-TMMC simulations require fixed values for the 
activity ^ 0], the particle-center accessible dimensions 
{hx, hy^ hz} that define the volume of the simulation cell 
V = hxhyhz, and the reciprocal temperature f3. For all 
simulations conducted here, we set f = 1, /? = 1, and we 
used the particle- wall interactions given by Eq. The 
values of hy x hz or hz are determined by the confining 
dimensions of the Id or 2d channels, respectively, and 
the remaining periodic dimcnsion(s) were chosen to sat- 
isfy V = 1000. Indistinguishable results were obtained 
for systems of size V = 500. 

The key quantities extracted from the GC-TMMC sim- 
ulations were the total particle number probability distri- 
bution n(A^), the excess configurational energy U°^{N), 
and the TV-specific spatial density distribution p{N,r), 
each evaluated over a range of particle numbers span- 
ning from = to iV = 984. Using basic arguments 
from statistical mechanics 0, 0| , one can relate these 
quantities to s'^^ for the inhomogeneous HS fluid: 

s°''(iV)/fcB =iV"^{ln[H(A^)/H(0)] -Mn^ + lnA^! 
-N\nN + PU''''{N)+ J p{N,r)\np{N,r)dr} (2) 

Given that V = 1000 is flxed, Eq. © provides s°''(/3/.) 
within the range < pft, < 0.984, where ph — N/V is the 
number density based on the particle-center accessible 
volume. 

First, we discuss our observations for the utility of ex- 
cess entropy s'^^ in predicting how confinement affects 
dynamics. Fig. ^ shows a parametric plot of D versus 
—s'^^ for the HS fiuid both in the bulk and confined to 
the 17 different 2d channels {hz = 2.5, 5, 10 with ew = 1, 
0.5, 0, -0.5, -1 and hz = 7.5, 15 with = 0). The data, 
which encompasses the dynamic behavior of the equilib- 
rium fluid from the dilute gas to the freezing transition, 
spans three decades in D. The collapse of the data onto 
a single master curve indicates that, to an excellent ap- 
proximation, the simple one-to-one correspondence be- 
tween D and s'^'^ for the bulk HS fluid also holds when 
the fluid is severely confined. Moreover, the quality of 
the data collapse is largely independent of either channel 
width (including "particle scale" channels with hz = 2.5) 
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FIG. 1: Self-diffusivity D versus the negative of excess en- 
tropy per particle — s^^/fce for the bulk HS fluid (solid curve) 
and for the HS fluid in 2d channels (symbols). The sym- 
bols correspond to hz — 2.5 (circle), 5 (square), 7.5 (plus), 
10 (triangle up) and 15 (cross). The color codes are =1 
(blue), 0.5 (cyan), (red), -0.5 (yellow), and -1 (green). Inset 
is the same plot with added points for the Id channels: 5x5, 
7.5 X 7.5, and 5 x 7.5 shown by red diamonds. 

and the sign or magnitude of the particle-wall interac- 
tion. Data for the HS fluid confined to the 3 rectanglular 
Id channels described above are superimposed in the in- 
set of Fig. n As can be seen, they also approximately 
collapse onto the bulk HS relationship between D and 
s™. In fact the diffusivities for approximately 50, 70, 
and 90% of the confined state points shown in Fig.^are 
within 3, 5, and 10 % respectively of the bulk value at 
the same excess entropy. 

Having established that s°^, an equilibrium thermo- 
dynamic property, can be used to predict how various 
confining environments impact the self-diffusivity of the 
HS system, it is natural to press forward and ask whether 
the same predictive information is contained in an even 
simpler structural measure: the average density. Here, 
one needs to be specific because there are two different 
definitions of average density that are commonly used 
to characterize inhomogenenous HS fluids {ph = N/V 
and p = N/V'). The former is based on the particle- 
center- accessible volume V — h^hyhz, while the latter 
is based on the total system volume; i.e., V = V for 
bulk systems, V = hxhy{hz + 1) for 2d channels, and 
V = hx(hy + ^){hz + 1) for Id channels. Schmidt and 
Lowen ji^ and, later, Zangi and Rice demonstrated 
that p is in fact the relevant density for the lateral (i.e., 
periodic) component(s) of the pressure, and thus p is also 
a natural independent variable for the free energy of the 
system. Below, we investigate whether p is also an accu- 
rate predictor for how conflnement effects HS diffusivity. 

Fig-Elshows D as a function of p for the bulk HS fluid 
as well as for the HS fluid confined to the 2d and Id chan- 
nels. The collapse of the data, while not perfect, demon- 
strates a very strong correlation between D and p that is 
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FIG. 2: Self-diffusivity D vs average density p — N/V' based 
on the the total system volume V' . Systems shown include 
the bulk HS fluid (solid curve) and the HS fluid confined to 
the 2d channels (symbols) discussed in the text. The inset 
includes the data from the HS fiuid confined to the Id chan- 
nels. Symbols for the main panel and the inset are identical 
to those given in Fig. Q 

nearly independent of the confining dimensions and the 
particle-wall interactions. Approximately 70 and 90% of 
the diffusivities for the confined state points shown in 
Fig. 12 are within 10 and 20%, respectively, of the bulk 
value at the same average density. This is another signif- 
icant result because, unlike s°^, p is intuitively simple to 
understand and trivial to determine in simulations (e.g., 
it is specified in microcanonical and canonical simula- 
tions). As is expected, the systems that exhibit the most 
noticeable deviation from bulk behavior in Fig.|21are the 
ones for which the fluid is under the most severe confine- 
ment, i.e., channels with dimensions comparable to the 
particle diameter. In these cases, it appears that spe- 
cific fiuid structuring (e.g., density enhancements in the 
channel corners) acts to only slightly reduce the diffusiv- 
ity relative to what would be predicted by the average 
density p. 

Given the approximate collapse of the data in Fig. |21 
it is natural to wonder whether it is simply the parti- 
cle structuring that is determining the HS dynamics. To 
test this idea further, we examine in Fig. 13 the local den- 
sity profiles p{z) for a HS fluid confined to 2d channels 
of width hz = 2.5 but with three different particle- wall 
interactions: = 1 (repulsive walls), e„ = (neutral 
walls), and tyj = —1 (attractive walls). All three systems 
exhibit the same average density p, and thus according to 
Fig. 121 display approximately the same self-diffusivity D 
as the bulk fluid. Clearly there are real and pronounced 
differences in the local structuring of the three confined 
fiuids, especially when compared to the uniform bulk ma- 
terial. These types of structural differences are usually 
the main focus of studies of inhomogenenous fluids by 
classical density functional theories. Interestingly, these 
pronounced structural differences only slightly alter both 
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FIG. 3: Density profiles p{z) for HS fluids conflned to a 2d 
channel of width hz — 2.5. Although each system has different 
particle- wall interaction strengths e^, they share a common 
average density p (dashed) and self-diffusivity D. 
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FIG. 4: Self-diffusivity D versus average density ph ~ N/V 
based on the the total particle-center-accessible volume V. 
Systems shown include the bulk HS fluid (solid curve) and 
the HS fluid conflned to the 2d channels between hard walls. 



s'^^ and D when considering fluids with the same average 
P- 

Finally, we examine what conclusions regarding den- 
sity follow if one instead chooses to plot D versus the al- 
ternative definition for average density, ph = N/V, based 
on the total particle-center-accessible volume V. In par- 
ticular, Fig. 21 compares the self-diffusivity for the HS 
fluid confined to 5 different 2d channels with hard walls 
(ew = 0) as a function of ph- Unlike when plotting versus 
p, there is no data collapse of D when plotting versus 
Ph- Thus, one might consider p a more natural indepen- 
dent variable than p/, . not only for thermodynamics of 
inhomogeneous fluids |20|,|2l|, but also for dynamics. 

To conclude, we have probed the structure, entropy, 
and diffusivity of the HS fluid confined to 2d and Id 
channels with a wide range of dimensions and particle- 
boundary interactions. Our main finding is that the re- 
lationships between diffusivity, excess entropy, and aver- 
age density for the bulk HS fiuid also remain valid, to 
within an excellent approximation, for the HS fiuid con- 
fined to particle-scale geometries. Since statistical me- 
chanical theories can provide accurate estimates for how 
confinement modifies the excess entropy and density, the 
robust connection between thermodynamics and dynam- 
ics reported here should have far-reaching implications 
for the prediction of dynamics in confined systems. We 
are currently testing whether similar connections hold (i) 
for the HS fluid in more general random environments 
(e.g., in quenched media, etc.), and (ii) for fluids with 
strong interparticle attractions that significantly affect 
local structuring. 
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